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Superconductivity is associated with spontaneously broken gauge symmetry. In some exotic 
superconductors the time-reversal symmetry is broken as well, accompanied with internal magnetic 
field. A time-reversal symmetry broken (TRSB) superconductivity without internal magnetic field 
involved can be induced by frustrated inter-component couplings, which becomes a realistic issue 
recently due to the discovery of iron-pnictide superconductors. Here we derive stability condition 
for this novel TRSB state using the Ginzburg-Landau (GL) theory. We find that there are multiple 
divergent length scales, and that this novel superconductivity cannot be categorized by the GL 
number into type I or type IL We reveal that the critical Josephson current of a constriction 
junction between two bulk superconductors of different chiralities is suppressed significantly from 
that for same chirality. This effect provides a unique way to verify experimentally this brand new 
superconductivity. 

PACS numbers: 74.50.-|-r, 74.25.Gz, 85.25.Cp 



Spontaneous symmetry breaking is one of the most 
fundamental concepts of modern physics, which plays a 
central role in particle physics, condensed matter physics, 
and even in our understanding into the birth of universe. 
Broken gauge symmetry is known as the hallmark of su- 
perconductivity 1, 2]. Simultaneous breaking of gauge 
symmetry and time-reversal symmetry has attracted con- 
siderable attentions Q, with the triplet superconductiv- 
ity in Sr2Ru04 as an icon Q- 

Possibility of time-reversal-symmetry-broken (TRSB) 
superconductivity in multi-component system was dis- 
cussed before in the context of conventional mechanism 
for exotic superconductivity jp . Interests in this novel 
phenomenon are renewed (a, |7| by the discovery of iron- 
based superconductors [sj. Previous works revealed that 
frustrated inter-component couplings induce intrinsically 
complex order parameters, and thus the TRSB super- 
conductivity ISMtI , primarily based on a special situation 
with all components equivalent (the isotropic system). 

It was also discussed that a junction structure between 
superconductors of two components and single compo- 
nent can exhibit a similar TRSB situation [9|. Possibil- 
ity of TRSB state was discussed in iron-based supercon- 
ductors with coexisting s- and d-wave order parameters 
[T^]. A novel dynamic mode was proposed for three- 
component, time-reversal symmetry reserved (TRSR) su- 



perconductivity . All these make the physics of multi- 
component superconductivity very rich. 

In the present work, we consider an anisotropic sys- 
tem with frustrated inter-component couplings. For sim- 
plicity, each component is considered as s-wave and well 
described by the BCS theory [ij. We adopt the GL ap- 
proach 12}, , powerful for discussions on stability of states 
and thus phase transition, symmetry breaking, as well as 
magnetic responses of superconductivity, which are the 
primary concerns of the present study. Benefitting from 
the simplicity of the GL theory, the condition for stable 
TRSB state are derived explicitly. It is also found that 




FIG. 1: TRSB states characterized by complex order param- 
eters and opposite chiralities of a three-component supercon- 
ductor in a Josephson junction of constriction structure. 



there are multiple divergent length scales, and that this 
novel superconductivity cannot be categorized by the GL 
number into type I and type II. We reveal that the critical 
Josephson current of a constriction junction between two 
bulk superconductors of different chiralities (see Fig. 1) is 
suppressed significantly from that for same chilarity. A 
standard Josephson junction measurement can provide 
smoking gun evidence for this brand new superconduc- 
tivity. 

The density of GL free-energy functional for a multi- 
component superconductor is: 

/ = E [%#.f + + 2^ l(f V - f A) 

' -E7jk(V'jV^:+c.c.) + gL(Vx A)2, 

j<k 

(1) 

where summations running over the number of compo- 
nents are understood, and all quantities are convention- 
ally defined (see 0, 0, HI, [H); aj{Tcj ) = specify the 
transition points for individual components. To be spe- 
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cific, we discuss here a system with three components, 
and extension of the following discussions to more compo- 
nents is straightforward. It is noticed that the physics ad- 
dressed below remains unchanged even when other pos- 
sible terms (up to the quartic order of order parameters) 
are included in Eq.(IT]). 

For 712723713 > 0, the system behaves basically as a 
single-band superconductivity close to the critical point 
of the composite system. Hereafter, we focus on the case 
7i27237i3 < 0, and treat a system with all jjk's negative, 
noticing that a simple gauge transformation in GL free 
energy ([T]) links it to other possible cases. The situation 
treated here may be realized in over-doped iron pnictides 

Around the critical point, the GL equations in absence 
of magnetic field can be linearized as 



(2) 



or in a vector form Q-^* = with coupling matrix Q. The 
critical point of the composite superconductivity Tc is 
given by the highest temperature where the determinant 
of Q becomes zero, i.e. 
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where temperature dependence of the quantities is under- 
stood (see Fig. 2). From Sylvester's criterion 1^, one has 



> and ajifc — 7jfc > at Tc, namely the critical point 
of the composite superconductivity is above any single 
components, and not below any of the two-component 
ones, in spite of negative couplings [T5| . 

If Eq.© has a single root at T = (see Fig. 2), 
or equivalently there are two independent vectors in the 
coupling matrix Q, the ratios among the order parame- 
ters given by the Cramer's rules for the components of 
the matrix Q should be real fT3|. The order parameters 
of three-component superconductivity that minimize the 
GL free energy ([T} can then be taken as real numbers, 
apart from a common phase factor, same as single- and 
two-component supernconductivity. 

Equation ^ has a doubly degenerated root [11] (see 
Fig. 2), or equivalently there is only one independent 
vector in Q, when 

0102 - 7i2 = 0> «ia3 - 7i3 = 0> 0203 - 7I3 = 0, (4) 

at Tc, namely Tc = Tci2 = Tci3 = Tc23. The single 
independent vector in the coupling matrix Q leaves the 
room for complex order parameters in Eq. ^ despite that 
all the parameters in Q are real. Equation Q is the first 
condition for a TRSB state specified by complex order 
parameters. 

It is easy to see that relations in Eq. , as well as the 
associated ones ajjki + Ijklji = 0, correspond to single 
zeros, since ^jk 7^ 0. One finds {akai — Jki)/i'^j^i~^ji) ~ 

(ofc ijkikih]i)l{aj + ijkijihki) = {ikihjif- 




FIG. 2: (color online). Schematic temperature dependence 
of determinant of the coupling matrix of the linearized GL 
equations (see Eqs.((2| and (O). The doubly degenerated root 
at Tc is a necessary condition for stable TRSB state, whereas 
the single-root case is associated with a conventional, TRSR 
state similar to single-component systems. 



The order parameters for T < Tc are given by 
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Putting -01 real as always possible, the imaginary parts 
in 02 and 03 should obey the relations 



(6) 



Therefore, for complex order parameters one has 
a362|02p + a263|03p ~ -a2a3 -I- 7I3 up to 0(i) with 
t = (Tc — T)/Tc. In the same way, one obtains two other 
similar relations, and thus 
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We then arrive at the following temperature dependence 
of order parameters: 
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up to 0(t). 

The single independent relation in Eq.([5]), for example, 
ai — 7i20'2/0'i — 7i3-03/0'i = 0, is then equivalent to 
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for T < Tc, where 02 1 (03i) is the phase difference be- 
tween (0^3) and ipi. It becomes clear that the condi- 
tion for a state of complex order parameters to be stable 
is equivalent to that of a triangle formed by three seg- 
ments: 
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FIG. 3: (color online). Phase diagram for a three-component 
superconductor with stable TRSB superconductivity at the 
central part and TRSR one at the corners. 




(a) 
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FIG. 4: (color online). Distortions in the complex order pa- 
rameters in TRSB state of the isotropic system (see text), 
(a) Definitions of distortions in order parameters, (b) and 
(c) Characteristic distortions of mode-I and model-II respec- 
tively, for m = m' . 
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for T < Tc. A phase diagram is displayed in Fig. 3. 

The relations @ and (|10p formulate the full condition 
for the stability of state with complex order parameters, 
i.e. the TRSB superconductivity. It is clear that the spe- 
cial case with isotropic parameters discussed previously 
[B-0] satisfies these conditions. 

A phase transition at a lower temperature Ttr < Tc 
from a TRSB state to a TRSR state is possible for 
appropriate temperature dependence of parameters (see 
Fig. 3), where interesting physics is expected. Although 
discussions in the GL scheme can be pushed forward, we 
notice that to treat the two transitions concretely one 
needs a more microscopic theory. 

Next we investigate the coherence length of the state 
of complex order parameters. In order to demonstrate 
the novelty of this state, we concentrate for a while on a 
system of isotropic parameters in Eq.(IT]) except mi = m 



and TO 2 



TO3 



the simplest, but non-trivial case 



allowing analytic treatment. In the bulk, the amplitude 
of the order parameter is V'o = IV'jl = -^Z— (a -I- 7)/&, 
and the phase difference between any of the two com- 
ponents is 27r/3 [5-7]. A local distortion in the first 
component '01 = (1 + ^i)V'o causes distortions in the 
other two components V'2 = (1 + <52)V'o exp[i(27r/3 -|- ^3)] 
and ip^ — [1 + ^2)0^0 exp[i(47r/3 — Ja)], with 5j's real, 
as displayed in Fig. 4a. In the one-dimensional case, 
the GL equation for the variation of ipi, for example, 
is oVji + &V'oV'i -7(^2 + V's) = (h^ /2rn)d'^il)i/ dx^ , which, 
with two other similar equations, yield 



The coherence length defined by 5j = exp(— -\/2a;/C) 
at large distance limit is thus determined as 
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There are two divergent solutions, and the corresponding 
characteristic modes are given by A2/A1 = — i?/[2(3 — 
2R)] and A3/A1 = -^3(2 - i?)/[2(3 - 2R)], where R 
is defined by the right-hand side of Eq. (IT^ . Mode-I as- 
sociated with the larger solution, thus giving the coher- 
ence length of the system, is specified by i? = 1, and 
^2Mi = -1/2 and A3/A1 = V^/2 when m = m' (see 
Fig. 4b), whereas mode-II by i? = 2, and A2/A1 — 1 
and A3/A1 — (see Fig. 4c), with the distortion vec- 
tors form an equilateral triangle in both cases. While 
mode-II at TO = to' is the conventional one associated 
merely with variation of amplitude (see Fig. 4c), known 
in single- and two-component cases, mode-I is the novel 
one in which variations of amplitude and phase are cou- 
pled (see Fig. 4b), specific to the TRSB state. 

Other quantities for the present superconductivity are 
available straightforwardly, such as the London penetra- 
tion depth A 
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the thermodynamic critical field Htc from Eq.([S 
1 

= 2(l + 7i27i3/ai723)^Xl"i/^J' 

j 

and the nucleation field (or Hc2) 
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/ = «i sin 93 + 12 sin(v3 — 202i) + «3 sin(93 — 203i), (16) 



It is clear that Htc 0o/2V27r^A, and iJn ^ 0o/27r^2. 
Consequently, we are led to conclude that this super- 
conductivity cannot be categorized by the GL number 
K = A/^ into type I or type II [l^ . 

Because of the existence of two divergent solutions of 
Ea.([T^. a vortex can exhibit different core sizes [12|] for 
the three components even close to T^. This may cause a 
long-range attractive and short-range repulsive interac- 
tion between two vortices and thus exotic vortex states, 
a possibility discussed previously in two-component case 

m 

Finally let us investigate the Josephson current be- 
tween two bulks linked by a short, narrow constriction in 
otherwise continuous superconducting material, as dis- 
played in Fig. 1. The length of the constriction is much 
shorter than the coherence length / <C ^. Suppose that 
the two bulks exhibit opposite chiralities by chance in 
a cooling process, with ifiji^ and "^jR the wave functions 
at the left and right bulk respectively, where |'0jL| — 
IV'jrI = V'jo, V'ir/V'il = e'"^, V2l/V'il = e'^^^ V2o/V'io, 
V'sl/V'il = e^'^^^Vso/V'io, ^2r/V'ir = e'^^— <A2i)^2o/V'io, 
and V'3r/V' 



with ij = eti^ip'jQ/mj. We arrive at the critical Josephson 
current 
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'03o/V'io- The order parame- 
ters on the bridge ipj{x) = fj(x)'ipjh are determined by 

ai/l + &lV'lol/l|Vl - 712/2'02l/'01L - 713/3V'3l/V'1L = 

{ti^/2mi)d'^fi/dx'^, and two other similar equations. 
Following the idea developed for conventional single- 
component superconductors [ip, [l^ ) one gets the linear 
solutions /i(x) = (1 - x/l) + e^'^x/l, f2{x) = (1 - x/l) + 
gi(2,r-202i+<p)^/; and/3(a;) = {l-x/l)+e'^^-^~^'>'''^+^^x/l, 
with a; = 0, Z at the left and right ends of the link, since 
{h'^/2mj)d'^fj/dx'^ = O(V'fo) cx C"^ for / < ^ 

The Josephson current between the two bulks is then 




J jik cos 2(j)jk, 

3<k 



(17) 



which equals zero in the isotropic case. On the other 
hand, when the two bulks share the same chirality, one 
has / = (ii -I- 12 +13) sinyj in contrast to Eq. ([T6l) . 

Therefore, because of the interference among the compo- 
nents, the critical Josephson current between two bulk 
superconductors of different chiralities is much smaller 
than the one for same chirality. 

In experiment based on the constriction structure 
shown schematically in Fig. 1, the two bulk superconduc- 
tors can show either the same chirality or the opposite 
ones by chance in repeated cooling processes. Then the 
standard measurement of critical Josephson current will 
give two sets of values, quite different from each other as 
discussed above, which gives a smoking gun evidence for 
this novel superconductivity. 

To finish we notice that although the GL approach is, 
in principle, justified only close to the critical point, the 
phenomena revealed in present are expected for the whole 
temperature regime. 
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